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INTRODUCTION 


The  recent  interest  in  the  near-millimeter-wave  region  of  the  electromag¬ 
netic  spectrum  has  created  a  need  for  attenuators,  phase  shifters,  isolators, 
and  other  devices  at  these  frequencies.  Attempts  to  scale  standard  ferrite 
devices  from  the  X  band  (9  GHz)  to  this  frequency  range  are  frustrated  by  the 
requirement  of  high  static  magnetic  fields.  At  X  band  the  magnetic  field 
required  for  ferromagnetic  resonance  is  approximately  3,000  gauss  (G*)^  while 
at  90  GHz  the  field  required  for  ferromagnetic  resonance  is  approximately 
30,000  G.  Such  high  fields  are  difficult  to  generate  for  use  in  practical, 
portable  devices.  One  method  of  overcoming  this  obstacle  is  to  use  ferrites 
with  large  internal  fields.  Then  the  resonance  can  be  achieved  by  applying  a 
modest  external  field. ^  However,  such  ferrite  materials  generally  have  a 
relatively  large  intrinsic  loss  and  this  leads  to  devices  that  have  a  large 
and  undesirable  insertion  loss.  We  wish  to  develop  a  new  technique  of  using 
conventional,  low-loss  ferrite  material  in  an  unconventional  manner  to  produce 
the  necessary  devices  at  near -mi Hi meter  wavelengths. 

In  this  report  we  investigate  the  possibility  of  building  a  low  applied 
magnetic  field  isolator  or  nonreciprocal  phase  shifter  at  90  GHz  or  at  higher 
frequencies  using  standard,  low- loss  ferrites.  The  system  studied  is  shown  in 
figure  1 .  It  consists  of  a  thin  ferrite  slab  backed  by  a  metal  plate.  The 
static  applied  magnetic  field  is  in  the  plane  of  the  ferrite  slab.  Solutions 
to  Maxwell's  equations,  including  the  equation  of  motion  of  magnetization,  are 
given.  These  equations  are  applied  to  the  particular  geometry,  and  it  is 
shown  that  under  appropriate  conditions  only  the  dominant  surface  wave  mode 
propagates.  The  dominant  surface  wave  is  the  transverse  electric  (TE)  surface 
mode  with  the  electric  vector  parallel  to  the  static  magnetic  field. 
Conditions  on  the  slab  thickness  are  found  such  that  only  the  dominant  TE  mode 
propagates.  As  the  static  magnetic  field  is  increased  from  zero,  the  phase  of 
the  forward  wave  (waves  propagating  along  the  positive  y-axis),  differs  from 
the  phase  of  the  backward  wave  (waves  propagating  along  the  negative  y-axis), 
and  this  can  be  used  as  a  nonreciprocal  phase  shifter.  If  the  static  magnetic 
field  is  increased  and  an  appropriate  slab  thickness  is  chosen,  one  of  the 
waves,  either  the  forward  or  reverse  wave,  is  cut  off,  and  this  can  be  used  as 
an  isolator.  Furthermore,  conditions  on  the  slab  thickness  and  static  mag¬ 
netic  field  can  be  found  such  that  neither  forward  nor  reverse  waves 
propagate.  The  field  strength  can  then  be  changed,  and  one  of  the  waves  will 
propagate.  Thus,  by  changing  the  applied  magnetic  field  suddenly,  we  have  a 
switching  device. 


^S.  Lax  and  K.  J.  Button,  Microwave  Ferrites  and  Ferrimagnetics,  McGraw 
Hill,  New  York  (1962),  [This  book  contains  a  large  number  of  practical 
applications  of  ferrites  to  microwave  devices.  Much  of  the  pioneer  work  on 
yttrium  iron  garnet  single  crystals  and  applications  of  ferrites  to  microwave 
devices  was  done  at  Harry  Diamond  Laboratories,  Reference  to  this  work  is 
contained  in  the  selected  bibliography  section,] 

^ff,  H,  Von  Aulock,  Handbook  of  Microwave  Ferrite  Materials,  Academic  Press, 
New  York  (1965),  p  451, 
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T  x  =  a 
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2.  THEORY 


In  general,  Maxwell's  equations 
(in  cgs  units)  when  the  fields  vary 
harmonically  as  e”^***^  are^ 


VxB  =  ikB 


VxH  =  -ikD  +  J 
c 


=  0 


7»D  =  4Tip 


k  =  u)/c 


I  For  our  purposes,  we  shall  assume  that 

IVIEJAL  ^  p  are  zero;  thus,  the  electric 

field,  B,  and  the  electric  displace- 
Figure  1 .  Coordinate  system  for  two  ment,  D,  are  related  by 
geometries  applicable  to  present 

analysis.  D  =  eS  .  (2) 

With  these  assumptions  and  equation  (2),  we  have 
7»B  =  0 


7xh  =  -ikeB  . 

The  configurations  we  have  chosen  are  shown  in  figure  1.  The  analysis  is 
equally  applicable  to  the  configuration  of  figure  la  or  figure  1b.  We  assume 
that  none  of  the  fields  vary  in  the  z  direction,  and  we  assume  a  TE  mode  in 
the  structure  (E^  =  E  =  0,  E^  *  0).  With  these  assumptions,  the  first  of 
equations  (1)  can  be  written  as 

ly  =z  = 


^E 

3x  z 


^J.  D,  Jackson,  Classical  Electrodynamics,  second  edition,  Wiley,  New  York 


In  appendix  A  we  obtain  the  B,  H  relationship 


Bx  =  +  iKHy 


and 


(A-13) 


By  =  WHy  -  ixH^  , 

When  these  relations  are  used  in  equation  (4)  we  get 
ly  =  ik(pH^  +  ixHy) 

and  ( 5 ) 

-  i<«x)  • 

The  results  given  in  equation  (5)  can  be  used  to  determine  the  magnetic  fields 
as 

«x  =  ly  ®  Ez)/[ik(li2  -  <2)] 

and  ( 6 ) 


Equation  (6)  determines  the  magnetic  field  components  once  the  electric  field, 
E^,  is  known.  To  determine  E^,  we  must  obtain  the  wave  equation  satisfied  by 
this  field  component. 

To  obtain  the  wave  equation  for  E^  we  take  the  curl  of  the  first  of 
equations  ( 1  ), 

VxVxE  =  ikVxB 


and  in  Cartesian  coordinates 

V(V.E)  -  v2e  =  ik(VxB)  . 


From  equation  (3)  and  since  E  has  only  a  z  component,  we  obtain 
=  Miy  Bx  -  h  ®y) 


where 


The  result  on  the  right  side  of  equation  (7)  can  be  reduced  by  using  equation 
{A-13)  and  V»B  =  0  to  obtain 


(8) 
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This  is  the  wave  equation  satisfied  by  for  the  region  occupied  by  the 
ferrite  slab.  To  obtain  the  corresponding  wave  equation  for  the  exterior 
region  (x  >  a),  let  k  =  0  and  e  =  p  =  1  in  equation  (8);  thus 

=  0,  X  >  a  .  (9) 

If  we  assume  that  we  have  a  wave  propagating  along  the  y-axis  so  that  all  the 
fields  vary  as  e^^^,  then  equations  (8),  (9),  and  (6)  can  be  written  as 


0  <  X  <  a  , 


X  >  a  , 

=  (quE  -t-  <  '  0  <  X  <  a  , 


H  =2-e 

X  k 


X  >  a 


(10) 


H..  = 


q<E)/[ik(u2 


-  .^2 


)] 


0  <  X  <  a  , 


and 


H..  = 


zl  dE 

ik  dx 


X  >  a  , 


where 


=  Ee 


iqy 


If  the  magnetic 


The  nonreciprocal  terra  appears  in  H^,  in  the  product  qic. 
field  is  reversed,  =  -<(8^]  (see  eq  (A-14)  and  (A-11),  app  A),  then  H 


returns  to  the  same  form  if  we  reverse  the  direction  of  propagation.  That  is. 


q<(Ho)  =  (^)'c(-H^) 


(11) 


However,  the  results  given  in  equations  (11)  and  (10)  are  not  sufficient  to 
assure  nonreciprocity.  To  achieve  a  truly  nonreciprocal  device  it  is  neces¬ 
sary  to  have  some  asymmetry  in  the  x  direction.  For  the  case  considered  here, 
the  asymmetry  is  supplied  by  the  conductor  at  x  =  0.  Asymmetry  could  also  be 
supplied  by  replacing  the  conductor  by  a  dielectric  slab.  The  solution  for  a 
dielectric  slab  is  considerably  more  complicated  than  for  the  conducting  plate 
considered  here,  but  the  problem  using  a  dielectric  slab  may  be  more  useful 
for  a  practical  device. 


We  wish  to  have  surface  waves  propagating  along  the  ferrite  slab,  so  we 
assume  solutions  for  equation  (10)  of  the  form 


0  <  X  <  a 


E  =  A  sin  (r^x)  , 


and  (12) 

E  =  ,  X  >  a  , 

where 

r?  =  k2e[u2  _  K^]/yi  -  q2 

and 


The  corresponding  magnetic  fields  are  obtained  using  equation  (10).  The 
boundary  condition  at  x  =  0  is  that  =  0,  and  the  solution  chosen,  in 

equation  (12),  assures  that  this  boundary  condition  is  satisfied.  The  bound¬ 
ary  conditions  at  the  surface  of  the  slab  (x  =  a)  are  that  E^  a  are 

continuous.  Using  equations  (12)  and  (10),  and  applying  these  lx>u.  ty  con¬ 
ditions,  we  obtain 

-(’o  =  (^("i  *  (1^) 

This  is  the  transcendental  equation  determining  the  propagation  cons  >  q,  if 
the  relationships  of  q  with  and  given  in  equation  (12)  are  usea. 

In  our  discussion  of  equation  (11)  we  pointed  out  that  the  nonreciprocity 
arises  due  to  the  product  qx.  This  is  also  true  in  equation  (13).  If  we 
express  P^  and  P^  in  terms  of  q,  using  equation  (12),  the  only  un)tnown  in 
equation  (13)  is  q.  All  the  quantities  in  equation  (13)  are  even  functions  of 
the  magnetic  field  (u,  P^,  and  P^),  except  that  k  is  an  odd  function  of  mag¬ 
netic  field  (<[-H^]  =  -k(h^)).  Thus,  if  we  find  the  solution  of  equation  (13) 
for  positive  (q^)  and  the  solution  of  negative  (q_)  assuming  q  positive 
in  both  cases,  then  Aq  =  q+  -  q_  would  be  the  differential  phase  shift.  The 
same  result  could  be  obtained  by  assuming  first  q  positive,  then  q  negative, 
but  the  magnetic  field  >  0  in  both  cases.  Because  of  the  complicated 
manner  in  which  the  magnetic  field  enters  in  equation  (13),  we  would  expect 
this  differential  phase  shift,  Aq,  to  vary  with  the  strength  of  the  static 
magnetic  field,  H  .  Thus,  by  varying  the  applied  magnetic  field,  H  ,  the 
differential  phase  shift  can  be  varied.* 


3.  COMPUTATION 

For  the  purposes  of  numerical  evaluation  it  is  convenient  to  introduce 
several  dimensionless  variables 


(fi.  Lax  and  K.  J.  Button,  Microwave  Ferrites  and  Ferrimagnetics ,  McGraw 
Hill,  New  York  (1962).  [Detailed  descriptions  of  various  nonreciprocal  phase 
shifters  are  given  in  the  later  chapters  of  this  work.] 


and 


(14) 


q  =  ko 

Using  these  variables  introduced  by  equation  (14)  in  equation  (13)  we  obtain 

-Yq  =  +  UYj^  cot  (kaYj^)]/(M2  -  ,  (13) 

and  equation  (12)  becomes 


and  (16) 


Y?  =  e[)i2  -  . 

The  solutions  of  equations  (15)  and  (16)  then  give  the  dimensionless 
quantities  Yj^»  and  a  for  any  particular  ferrite.  Because  of  the 

complicated  nature  of  these  equations,  it  is  extremely  difficult  to  obtain  a 
solution,  particularly  when  magnetic  losses  and  dielectric  losses  are 
considered.  We  shall  restrict  ourselves  to  the  lossless  case.  Losses  in  a 
dielectric  slab  (H^  =  0  in  eq  (15))  are  considered  in  appendix  B. 

The  lossless  case  for  the  magnetic  part  is  given  by  letting  a  =  0  in 
equations  (A-1 1 )  and  (A-14).  If  we  assume  that  the  dielectric  constant  is 
real  (to  consider  losses  in  the  dielectric,  we  assume  e  =  e'  +  ie",  with  e"  > 
0),  then  the  roots  of  equation  (16)  are  all  real.  A  useful  result  can  be 
obtained  from  equation  (15)  if  we  investigate  the  li-^it  Yq  0^  (that  is, 
Fq  0'*’),  This  situation  corresponds  to  the  limit  where  the  surface  wave  on 
the  ferrite  is  minimally  bound.  If  we  let  =  0  in  equation  (15)  we  have 

<+PYi^  cot  (jtay^)  =  0  ,  (17) 

and  for  Yq  =  0  equation  (16)  yields 

0=1  , 

(18) 

Y,  -  le(„2  -  k2)/„  -  1]’/^  , 


and  since  equation  (17)  is  even  in  y^,  we  need  only  consider  solutions  for 
which  y^  >  0.  The  solution  of  equation  (17)  for  Kay^  is  given  by 


(19) 


where  p  =  0,  1,  2,  .  .  .  (Actually  p  can  take  on  negative  values,  but  for  our 
purposes  p  >  0  is  sufficient.)  The  value  of  p  in  equation  (19)  is  the 
smallest  integer  that  makes  the  right  side  positive  (larger  positive  values  of 


10 


p  give  higher  modes).  Since  Yj^  is  given  by  equation  (18),  equation  (19)  gives 
the  thickness,  a,  of  the  ferrite  slab  for  a  particular  value  of  magnetic 


field,  Hq.  Thus  we  have 


and  if  we  let  +  -H^  we  have 


since  =  -<(Hq) 


.  -  (-  *  p+’]  ' 

we  have 

a_  =  ^  [tan"’  +  p_tt]  , 


For  the  lossless  case  (a  =  0)  using  equations  (A-14)  and  (A-11),  we  can 
write  the  quantities  entering  equations  (21)  and  (18)  as 


2  2 
V  -  <  = 


(“,  -  “)(«,  t  b) 


“  ("l  -  “2) 


u  .  1  + 


-  h: 


where 


B  =  H  +  4ttM  , 

H2  =  /bh  , 

=  OJ/Y  . 

The  analytical  form  corresponding  to  a  Brillouin  function**  for  spin  1/2  parti¬ 
cles 

M(H)  =  M  tanh  (sH)  ,  (22) 

o 

is  assumed  with  s  chosen  by  assuming  M(H)  to  be  90  percent  of  the  saturation 
magnetization,  at  a  field  of  500  Oe.* 

The  slab  thickness,  a^,  for  magnetic  fields  in  the  positive  z  direction 
given  in  equation  (20)  with  the  smallest  value  of  p_|^  gives  the  lowest  mode 
cutoff.  That  is,  any  slab  thinner  than  a^  will  be  cut  off.  The  higher  modes 

**0.  Kittel,  Introduction  to  Solid  State  Physics,  fourth  edition,  Wiley,  New 
York  (1971),  p  505. 

*(0e)  =  (A/m)  X  79.577 


are  given  by  larger  values  of  p_j_  in  equation  (20).  The  slab  thicKness  given 
in  equation  (21)  is  for  the  nvagnetic  field  in  the  negative  z  direction.  The 
resulting  thicknesses,  a_j_  and  a_,  f'  the  lowest  mode  for  yttrium  iron  garnet 
(YIG)  are  shown  in  figure  2.  The  frequency  chosen  is  90  GHz  and  the 
properties  of  YIG  are  from  Von  Aulock.^  From  figure  2  we  see  that  at  = 
3000  Oe,  for  a  slab  thicker  than  283  pm,  radiation  propagates  in  both  direc¬ 
tions.  With  =  3000  Oe  and  for  a  slab  between  283  and  276  pm,  only  waves 
propagating  in  the  positive  y  direction  are  allowed  and  we  have  an  isolator. 
For  slabs  less  than  276  pm  and  a  field  of  3000  Oe,  no  wave  propagates.  There¬ 
fore,  curves  as  given  by  figure  2  can  be  used  to  design  isolators,  switches, 
and  nonreciprocal  phase  shifters  for  ferrite  materials,  provided  the  material 
has  low  dielectric  and  magnetic  losses.  Even  in  cases  where  these  losses  are 
not  negligible,  data  as  presented  in  figure  2  are  helpful  for  design  purposes, 
provided  operating  points  are  chosen  as  far  as  possible  from  the  cutoff 
points.  Figure  3  shows  the  slab  thicknesses  a^  and  a_  for  a  nickel  zinc 
ferrite  at  90  GHz.  The  slab  thicknesses  are  reduced  because  of  the  larger 
dielectric  constant  (e  =  12.5  as  compared  to  e  =  10  for  YIG).  The  differen¬ 
tial  thickness  given  by 

Aa  =  a_  -  a^  (23) 

is  increased  because  of  the  larger  saturation  magnetization  (4ttM^  =  1750  G  for 
YIG,  4ttM^  =  5000  G  for  nickel  zinc  ferrite).  In  both  figures  2  and  3,  the  low 
applied  magnetic  field  values  (H^  <  500  Oe  in  fig.  2  and  H^  <  1000  Oe  in 
fig.  3)  should  not  be  given  too  much  consideration. 


Ho(0  ) 

Figure  2.  Slab  thickness  as  a  function  of  magnetic  field  for 
with  reported  values  at  X-band  frequencies. 
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Figure  3.  Slab  thickness  as  a  function  of  magnetic  field  for  nickel-zinc 
ferrite. 

To  illustrate  more  clearly  the  effect  of  changing  only  one  of  the  vari¬ 
ables  of  a  material,  such  as  the  dielectric  constant  or  saturation  magnet¬ 
ization,  several  hypothetical  cases  are  considered  using  the  values  for  YIG  as 
a  starting  point.  In  figure  4  the  dielectric  constant  of  the  ferrite  is 
increased  to  12  from  the  original  value  (10)  used  in  figure  2.  The  results 
show  that  thinner  slabs  are  required  for  similar  behavior,  but  that  the  dif¬ 
ferential  thickness.  La,  remains  essentially  unchanged  from  its  original  value 
in  figure  2.  Increasing  the  dielectric  constant  to  14  has  essentially  the 
same  effect,  as  can  be  seen  in  figure  5.  With  the  dielectric  constant  at  10 
but  with  the  saturation  magnetization  increased  from  1750  (in  fig.  2)  to  2000 
G,  the  calculations  are  repeated  and  the  results  are  shown  in  figure  6. 
Comparing  these  results  with  those  in  figure  2,  we  see  that  there  is  only  a 
slight  change  in  the  differential  thickness,  Aa.  Since  the  saturation  magnet¬ 
ization  changes  with  temperature,^  the  results  shown  in  figure  6  indicate  that 
an  isolator  would  not  be  very  sensitive  to  changes  in  the  ambient  temperature. 
Figures  7  and  8  have  larger  departures  of  the  saturation  magnetization  from 
the  value  of  1750  G  used  in  figure  2.  They  show  that  the  differential  thick¬ 
ness  increases  with  saturation  magnetization,  but  the  slab  thickness  is 
essentially  constant  for  changes  in  4TrM.  Therefore,  the  results  given  in 
figures  2  through  8  would  indicate  that,  to  a  good  approximation,  we  can 
assume  that  increasing  the  dielectric  constant  decreases  the  required  slab 
thickness  and  increasing  the  saturation  magnetization  increases  the  differ¬ 
ential  thickness.  These  two  effects  are  independent. 
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Figure  4.  Slab  thickness  as 
a  function  of  magnetic  field 
for  a  hypothetical  YIG  sample 
(dielectric  constant  assumed 
to  be  12). 
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Figure  6.  Slab  thickness  as 
a  function  of  magnetic  field 
for  a  hypothetical  YIG  sample 
(saturation  magnetization 
assumed  to  be  2000  G). 
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Figure  7,  Slab  thickness  as 
a  function  of  magnetic  field 
for  a  hypothetical  YIG  sample 
(saturation  magnetization 
assumed  to  be  2500  G). 


Figure  8.  Slab  thickness  as 
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a  function  of  magnetic  field  (2 
for  a  hypothetical  YIG  sample  g 
(saturization  magnetization  ^  279.2 
assumed  to  be  3000  G).  3 
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Next,  the  frequency  sensitivity  of  the  solutions  to  equations  (20)  and 
(21)  are  investigated.  For  the  purpose  of  this  investigation  YIG  was  chosen 
with  e  =  10  and  4irM^  =  1750  G  as  in  figure  2.  Figure  9  shows  the  variation  of 
slab  thickness,  a^,  and  differential  thickness,  Aa,  as  a  function  of  magnetic 
field  for  three  different  frequencies  (80,  90,  and  100  GHz).  In  figure  10  the 
thickness,  a.,  is  plotted  versus  frequency  from  80  to  100  GHz  for  a  magnetic 
field  Hq  =  3000  Oe.  The  slab  thickness  for  both  a^  and  a_  decreases  with 
increasing  frequency.  Thus,  if  the  dielectric  constant  decreases  with  fre¬ 
quency  in  the  correct  manner,  the  slab  thickness  would  remain  constant. 
Unfortunately,  ferrites  with  dispersion  in  any  region  of  interest  have  a 
prohibitive  amount  of  loss  in  that  region. 


4.  CONCLUSION 


We  have  found  the  conditions  under  which  a  surface  wave  propagates  on  a 
ferrite  slab  backed  by  a  metal  conductor.  The  only  region  investigated  was 
that  with  the  static  magnetic  field  much  less  than  the  magnetic  field  required 
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Figure  9.  -Slab  thickness  as  a 
function  of  magnetic  field  for 
YIG  for  three  frequencies  with 
reported  values  of  e  and  4tiM 
at  X  band. 
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Figure  10.  Slab  thickness  as 
a  function  of  frequency  for 
YIG  with  reported  values  of  e 
and  4TrM  at  k  band. 
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for  ferromagnetic  resonance  in  the  material.  The  results  reported  here  are 
for  the  dominant  TE  mode  of  the  system.  Magnetic  and  dielectric  losses  were 
assumed  to  be  negligible  in  this  analysis.  It  is  expected  that  these  results 
will  be  useful  in  the  design  of  a  practical  device,  provided  they  are  used  in 
a  region  some  distance  from  an  abrupt  cutoff.  The  abrupt  cutoff  points  should 
be  softened  by  the  introduction  of  losses  in  a  manner  similar  to  those  used  in 
waveguides  when  the  conductivity  of  the  walls  is  included. 

The  results  of  the  analysis  are  presented  in  the  form  of  graphs  that 
border  the  three  regions:  both  forward  and  reverse  propagation  but  with 
different  phase  velocities,  only  forward  propagation,  and  no  propagation.  At 
any  given  frequency  these  regions  are  outlined  by  slab  thickness  as  a  function 
of  applied  magnetic  field.  The  forward  and  reverse  propagation  in  the  first 
region  suggests  its  use  as  a  phase  shifter,  and  the  second  region  can  be  used 
for  isolation.  The  third  region  can  be  used  for  switching  if  the  applied 
magnetic  field  is  suddenly  changed  from  this  region  to  region  two. 
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Losses  were  considered  only  for  the  particular  case  where  the  applied 
magnetic  field  was  zero  and  magnetic  losses  were  ignored.  Small  dielectric 
losses  were  introduced  and  the  complex  propagation  constant,  q  =  q'  +  iq", 
along  the  surface  was  calculated.  It  was  found  that  for  slab  thicknesses 
restricted  so  that  only  the  dominant  mode  propagates,  the  value  of  q"  was 
smallest  for  the  minimum  slab  thickness.  A  maximum  occurred  in  q"  near  the 
middle  range  of  slab  thickness,  with  q“  gradually  decreasing  as  the  thickness 
of  the  slab  increases. 

The  predicted  slab  thickness  for  YIG  falls  in  the  range  of  thick  films 
which  can  be  fabricated,  and  perhaps  it  could  be  deposited  on  the  conductor  by 
a  method  similar  to  liquid  phase  epitaxy.  If  single-crystal  films  are  de¬ 
sired,  the  film  probably  would  have  to  be  deposited  on  a  nonmagnetic  single¬ 
crystal  substrate,  such  as  yttrium  gallium  garnet,  where  the  crystal  cell  size 
matches  the  magnetic  film  cell  size.  In  such  cases,  the  analysis  of  the 
system  would  have  to  be  altered  considerably  from  the  one  presented  here. 

In  further  analysis,  magnetic  as  well  as  dielectric  losses  in  the  ferrite 
slab  should  be  considered.  The  inclusion  of  these  losses  would  give  the 
insertion  loss  for  both  the  phase  shifter  and  isolation  mode  of  operation. 
Further,  for  an  isolator,  the  degree  of  isolation  could  then  be  estimated  much 
more  realistically. 
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The  equation  of  motion  of  magnetization,  as  given  by  Landau  and 
Lifshitz, ^  is 


ft'  =y(M'xH')  -  [m'x(||'xH' )  ]  ,  (A-1  ) 

where  y  is  the  gyromagnetic  ratio  and  a  is  the  damping  constant.  The  other 
constants  in  the  second  term  on  the  right  in  equation  (A-1 )  are  added  for 
convenience  so  that  a  is  dimensionless.  In  most  of  the  materials  that  we 
consider,  the  damping  is  small  and  it  is  convenient  to  use  Gilbert's 
approximation.^  That  is,  in  the  second  term  on  the  right  we  assume  that 

M'xH'  =  A'  .  (A-2) 

Y 

The  equation  (A-1 )  becomes 


A'  =  y(M'xH')  -  j^rj  (M'xA')  .  (A-3) 

The  result  given  in  equation  (A-3)  agrees  with  the  result  of  equation  (A-1 ), 
if  terms  of  order  and  higher  are  ignored.  At  the  low  power  levels  that  we 
considered,  it  is  sufficient  to  linearize  equation  (A-3)  by  letting 


and 

with  the  conditions 
and 


U' 

=  ez»o 
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A 

In  this  equation  is  a  unit  vector  along  the  z-axis.  In  equation  (A-4),  we 
assume  the  static  magnetic  field,  and  the  static  magnetization,  M^,  to  be 
in  the  z  direction.  The  assumption  of  small  amplitude  signals  also  leads  to 
M"  and  H"  having  only  components  transverse  to  z. 


If  we  use  equation  (A-4)  in  equation  (A-3)  we  have 


A"  =  YfA"><e  H  +  e  M  xH"]  -  fe  M  xA"l  , 

^  ZO  ZO  '  M''ZO  ' 


(A-5) 


where  we  assume  that  |h|  =  M^.  We  assume  that  all  the  time-varying  quantities 
have  a  harmonic  time  dependence,  and  we  let 


^L.  Landau  and  L.  Lifshitz,  On  the  Theory  of  the  Dispersion  of  Magnetic 
Permeability  in  Ferromagnct ic  Bodies,  Physik  Zeitschrift  Sowjetunion,  ^ 
(1935),  153. 

^T.  A.  Gilbert,  Armour  Research  Foundation,  Rept.  No.  11  (January  25,  1955). 
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M"  = 

and  (A-6) 

H"  =  , 

with  M  and  H  time  independent  but  still  spatially  dependent.  Using  equation 
(A-6)  in  equation  (A-5)  and  expressing  the  results  in  Cartesian  components 
gives 

-iojM^  =  Y(MyHo  ■  '^o“y^  ■ 

and  (A-7) 

-ioMy  =  -  M^H^)  +  . 

To  obtain  the  relation  of  M  to  H,  it  is  convenient  to  use  the  circular 
polarized  quantities, 


M 


± 


«x  * 


and 


(A-8) 


H  =  H  ±  iH 

±  X  y 

Substituting  equation  (A-8)  into  equation  (A-7)  we  obtain 

+YM  H 

M  =  - -  .  (A-9) 

t  oj(1  ±  ia)  +  YHq 

In  the  cgs  system  the  relationship  of  B,  H,  and  M  is 

B  =  H  +  4iiM  (A-10) 


so  that  we  can  write 


^  = 


From  equations  (A-10)  and  (A-9) 
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± 


Y4nM 

_  _  o _ 

^  (*)(  1  ±  ia)  +  YH 

o 


(A-1 1 ) 


The  properties  of  ferromagnetic  materials  are  characterized  by  4itM^  and  a. 
However,  the  damping  parameter  a  is  replaced  frequently  by  the  linewidth, 
measured  at  constant  frequency  as  a  function  of  H^.  The  parameter  a  is  then 
related  to  the  linewidth,  AH,  by 


a 


(A-12) 
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where  AH  is  half  of  the  full  linewidth  at  half  maximum  of  the  absorption. 
Sometimes  the  full  Linewidth  at  half  maximum  of  the  absorption  is  reported; 
thus,  equation  (A-12)  must  be  used  with  caution. 

The  circular  polarization  form  of  the  magnetic  permeability  is 
convenient  for  electromagnetic  boundary  value  problems  when  the  ferrite 
material  has  circular  symmetry  about  the  static  magnetic  field,  H^.  However, 
for  our  problem,  the  permeability  tensor  relation  of  B  and  H  is  more 
convenient.  This  relation  is  usually  written  as 

and  (A-13) 

By  =  hHy  -  i<H^  . 


From  equation  (A-11)  we  find 


U  =  (u  +  M  )/2 

and  ■  (A-14) 

<  =  (y+  -  . 


For  the  case  where  a  is  small  and  the  static  magnetic  field  is  far  below 
resonance  (YB^  <<  oj),  we  have 


and 
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Uurinij  the  course  of  this  investigation  we  consider  only  Lossless 
material.  The  reason  for  this  assumption  is  the  difficulty  of  developiny  a 
(jeneral  technique  for  examining  the  roots  of  the  resulting  transcendental 
equation  which  determine  the  modes  of  propagation  (an  example  is  (13)  in  the 
main  text). 

In  this  appendix  we  examine  the  consequence  of  assuming  that  the  mag¬ 
netic  field  is  zero  and  the  only  source  of  losses  is  the  dielectric.  Further, 
if  we  assume  the  losses  are  small,  the  resulting  equations  for  determining  the 
propagation  constant  can  be  solved  by  a  perturbation  technique.  We  assume 
that  the  dielectric  constant  is  given  by 


e  =  e '  +  ie"  (B-1  ) 

and  further  assume  that  e"  <<  e'.  If  we  let  H  =  0  ( u  =  1 ,  <  =  0)  in  equation 
(13)  in  the  main  text,  we  obtain 

cot  =  -Yq  (B-2) 

with 

Y?  +  Y^  =  £  ~  1 
'i  'o 

and 


Now  if  we  assume 


and 


+  lYV 


Y  =  Y'  +  iY” 
o  o  o 


(B-3) 


o  =  o'  +  io" 


with  all  the  imaginary  parts  small  compared  to  the  real  parts,  then 

2 


and 


cot  (kaYj^)  =  cot  (kaY^j  -  ikaY^[l  +  cot^  [kaY|)] 

i(e"  -  2y^yV) 


Y'  +  iY"  =  Y'  + 
o  o  o 


(B-4) 


(B-5) 


Thus  the  real  part  of  equation  (B-2)  satisfies 

Y;  cot  (kaY^)  =  -Y^  ,  (B-6) 

so  that  the  real  part  of  Yi  satisfies  the  same  equation  as  equation  (B-2) 
when  Yj^  is  replaced  by  y^.  The  imaginary  parts  of  equations  (B-4)  and  (B-5) 
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satisfy 

kaY['[l  +  cot^  (<aYi)]  =  (e"  -  2y^Y[')/Yo  ' 

which  can  be  further  reduced  by  using  equation  (B-6)  and 


to  give 


(B-8) 


The  remainder  of  the  imaginary  parts  of  the  quantities  of  equation  (B-3)  can 
be  obtained  by  using  the  result  given  in  equations  (B-8)  and  (B-2)  to  give 


y;;  =  (eV2  -  Y^Y|‘)/Yi 


and 


(B-9) 


a"  =  YoY^o'  • 

Thus,  once  the  solutions  of  equation  (B-6)  are  determined,  the  quantities 
given  in  equations  (B-8)  and  (B-9)  can  be  determined. 


The  solution  for  equation  (B-6)  has  been  obtained  for  e’  =  10  and  e"  = 
0.01  for  the  dominant  mode.  The  region  of  interest  of  thickness  a  in  equation 
(B-6)  is  easily  shown  to  be 


(B-10) 


where  for  smaller  thickness  than  the  lower  limit  no  surface  wave  propagates, 
while  for  larger  thickness  than  the  upper  limit  both  the  dominant  mode  and 
higher  modes  propagate.  The  solution  for  equation  (B-6)  for  f  =  90  GHz,  e'  = 
10,  and  e"  =  0.01  for  values  of  a  in  the  interval  is  given  by  equation 
(B-10).  The  results  for  q"  (q"  =  ko")  are  shown  in  figure  B-1 .  Since  all  the 
yield  quantities  vary  as  e^^^,  a”  gives  the  decay  of  the  fields  in  the  di¬ 
rection  of  propagation  as  e“*^®"|y|.  The  results  show  that  for  thin  slabs  (a 
near  the  lower  limit  in  eq  (B-10))  the  attenuation  is  small.  As  the  slab 
thickness  increases,  the  attenuation  increases,  reaching  a  maximum,  and  then 
decreases  as  a  approaches  the  upper  limit  in  equation  (B-10).  Also  shown  is 
the  same  calculation  with  e"  =  0.02  and  all  other  quantities  unchanged. 

During  the  course  of  the  computation  all  the  imaginary  parts  in  equation  (B-3) 
were  monitored  and  found  to  be  small  compared  to  the  real  part  of  the  cor¬ 
responding  quantities.  Thus  the  results  are  valid  for  the  entire  range  of  a 
given  in  equation  (B-10). 
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Figure  B-1 .  Imaginary  part  of  propagation  constant,  q  =  q'  +  iq",  with  e 
10  and  e"  =  0.1,  or  at  90  GHz.  Slab  thickness  restricted  so  that  only 
dominant  mode  propagates. 
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